Estimation of covariance matrix from geochemical data with observations below detection limits     

Chang-Jo F. Chung 《Mathematical Geology》1993,25(7):851-865

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.  相似文献   

硅质碎屑砂岩中胶结物来源的新认识     

杨香华 《地质科技情报》1993,12(4):34-39

分析了浅埋藏砂岩的成岩环境,并从岩石学、矿物学、同位素及质量平衡估算方面论证硅质碎屑砂岩中的胶结物主要来源于砂岩浅埋藏期间的地下水,进而分析浅埋藏带内地下水的流动机制及早期胶结作用对晚期成岩作用的影响。  相似文献   

施工顺序对软土地坪复合地基的影响探讨     

魏海涛  蔡智 《岩土工程技术》2021,(1)

场地回填土与搅拌桩施工顺序将直接影响软土地坪复合地基的处理效果,其中沉降控制是影响地坪安全运行的关键。以近海软土地区工业厂房水泥土搅拌桩复合地基工程为例,探讨了复合地基水泥土搅拌桩、回填土两种不同施工顺序的影响因素及相应处理效果;通过数值模拟试验,计算分析了不同填土厚度下两种工况的地基固结总沉降、施工沉降、工后沉降。结果表明,采用“先土后桩”的施工顺序的施工期沉降大于“先桩后土”,但工后沉降大大减小,提高了软土地坪地基处理效果,有利于处理后地坪的安全运行,验证了前述分析结论,为合理设计方案的确定提供了重要的参考依据。  相似文献   

Estimation for partially observed Markov processes     

M. E. Thompson  T. N. Kaseke 《Stochastic Hydrology and Hydraulics》1995,9(1):33-47

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.  相似文献   

湖南骑田岭花岗岩及其暗色微粒包体的地球化学与壳幔岩浆的混合作用   总被引：12，自引：2，他引：12  

付建明  谢才富  彭松柏  杨哓君  梅玉萍 《地球学报》2006,27(6):557-569

骑田岭岩体是南岭地区燕山早期具幔源组分贡献的花岗岩典型代表。其主体岩性为角闪黑云二长花岗岩和(角闪石)黑云正长花岗岩,其中发育暗色微粒包体和由暗色矿物组成的团块或条带。暗色微粒包体具有岩浆混合的大部分岩相学证据。如包体的浑圆状外形、塑性形变、冷凝边、斜长石An的“双峰式”分布、似环斑长石、针状磷灰石等。包体属于准铝质(A/KNC=0.72～1.00,平均0.85)钾玄岩系列岩石,寄主岩石为准铝或弱过铝质(A/KNC=0.89～1.06,平均0.97)高钾钙碱性系列的岩石。二者在主量和微量元素上表现出岩浆混合成因的演化趋势。包体与寄主岩石的同位素组成具趋同性,它们的ISr和εNd(t)值分别为0.71041～0.71263、-6.9～-5.3和0.70854～0.71416、-9.2～-5.1,均表现出壳幔混源花岗岩类岩石的特点。包体K-Ar年龄为152Ma,与其寄主岩石的形成年龄(155～161Ma)接近,显示岩浆混合作用发生的时间大致为晚侏罗世早期。对包体及其寄主岩石产出的构造背景和地球化学特征的综合分析表明,该岩体中的暗色微粒包体是在伸展作用的大地构造背景下,上涌的幔源基性岩浆及其诱发的长英质酸性岩浆混合作用的产物。  相似文献   

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data     

G. Grünthal  R. Wahlström 《Natural Hazards》2006,38(1-2):159-176

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.  相似文献   

Re-evaluation of focal depths and source mechanisms of selected earthquakes in the Afar depression   总被引：1，自引：0，他引：1  

L. Hagos  H. Shomali  R. Roberts 《Geophysical Journal International》2006,167(1):297-308

We present a stepwise inversion procedure to assess the focal depth and model earthquake source complexity of seven moderate-sized earthquakes  (6.2 > M _w > 5.1)  that occurred in the Afar depression and the surrounding region. The Afar depression is a region of highly extended and intruded lithosphere, and zones of incipient seafloor spreading. A time-domain inversion of full moment tensor was performed to model direct P and SH waves of teleseismic data. Waveform inversion of the selected events estimated focal depths in the range of 17–22 km, deeper than previously published results. This suggests that the brittle–ductile transition zone beneath parts of the Afar depression extends more than 22 km. The effect of near-source velocity structure on the moment tensor elements was also investigated and was found to respond little to the models considered. Synthetic tests indicate that the size of the estimated, non-physical, non-isotropic source component is rather sensitive to incorrect depth estimation. The dominant double couple part of the moment tensor solutions for most of the events indicates that their occurrence is mainly due to shearing. Parameters associated with source directivity (rupture velocity and azimuth) were also investigated. Re-evaluation of the analysed events shows predominantly normal faulting consistent with the relative plate motions in the region.  相似文献   

A chronological assessment of Lake Okeechobee (Florida) sediments using multiple dating markers     

Shawn P. Schottler  Daniel R. Engstrom 《Journal of Paleolimnology》2006,36(1):19-36

Reconstructing recent limnological history often relies on lead-210 dating to accurately ascribe a chronology to a sediment profile. In Lake Okeechobee, Florida, a large, shallow subtropical lake that may experience severe mixing, multiple dating methods are required to confirm that conformable sedimentation has been preserved and that the assumptions of the ²¹⁰Pb method are satisfied. This study uses stratigraphic profiles of heavy metals, ¹³⁷Cs, PCBs and pollen as independent dating markers to validate the sediment chronology as determined by ²¹⁰Pb for three cores from the central mud zone of the lake. Unsupported ²¹⁰Pb and most dating markers show distinct concentration/depth profiles, suggesting that the sediments have not been severely mixed for at least the last 75 years. Onset and maximum activity of the radioisotope ¹³⁷Cs in the cores coincides with the ²¹⁰Pb-dated interval of 1945–1970. This agrees well with the known timing of atmospheric deposition of ¹³⁷Cs that resulted from above-ground nuclear testing during late 1940s until 1963. Sediment core profiles of atmospherically deposited metals such as Zn and Pb, which reflect regional increases during industrialization and decreases after regulation in the 1970s, exhibit expected concentration increases and peaks coinciding within 5–15 years of the predicted ²¹⁰Pb dates. Uranium, a contaminant in some phosphate fertilizers, shows large concentration increases at core depths dated to be about 1950 by ²¹⁰Pb, matching the intensification of agriculture after WWII. PCBs, which are expected to peak in the early 1970s, were measured in one core, and the observed peak corresponds to a ²¹⁰Pb date of about 1960. Pollen makers were unable to verify specific events, but increases in disturbance taxa and declines in native types correspond generally with the expected dates assigned by ²¹⁰Pb dating. Conformity between the ²¹⁰Pb defined dates and independent markers of < ±15 years confirm that Lake Okeechobee sediments do preserve a sequential and reliable stratigraphic history of the lake, useful for reconstructing past limnological conditions.  相似文献   

Tackling social exclusion in European neighbourhoods: experiences and lessons from the NEHOM project     

Anneli Kährik 《GeoJournal》2006,67(1):9-25

Disadvantaged urban neighbourhoods with their complex social, economic and physical situations have increasingly become the focus of governmental policies in the EU. The emphasis has shifted from improving physical qualities to building opportunities for social groups in these areas. While physical upgrade is easier to achieve, at least in the short run, initiatives often fail in outcomes on communities and local participation. This paper aims to improve understanding on how to tackle area-based social exclusion and create socially coherent communities. Building on the results of the NEHOM project, the focus is on examining how various interventions have impacted local communities and their opportunities for better inclusion. A comparison between neighbourhoods shows that mixing social groups to a certain degree is likely to be a necessary step towards opening long-term opportunities for a disadvantaged neighbourhood. The planned initiatives should at the same time directly benefit local communities, and the involvement of local groups in the planning and implementation process helps to guarantee that the overall changes are accepted by the local community. However, the proposed guidelines do not provide a universal recipe for a successful renewal suiting each neighbourhood—the contextual matters should always be a starting point when planning an initiative.  相似文献   

Evidence of reworking and resuspension of carbonates during last glacial maximum and early deglacial period along the southwest coast of India     

Kamlesh Verma  M Sudhakar 《Journal of Earth System Science》2006,115(6):695-702

A gravity core collected from the upper slope of southwest of Quilon at a water depth of 776 m (Lat: 8°12′263″N, Long: 76°28′281″E) was analysed for texture (carbonate free), calcium carbonate and organic carbon. Variation in silicic fraction seems to be controlled by silt, i.e., enrichment from 15 ka BP to 10 ka BP and then constant in Holocene. Below 15 ka BP, the silicic fraction gets depleted compared to the Holocene section with a minimum around 21 ka BP. Clay content remains nearly constant except in the Holocene where it shows an enrichment. Carbonate content of less than 63 micron when computed by subtracting coarse fraction content from the total carbonate suggests that the total carbonates are mainly concentrated in the finer fraction. All these carbonate phases show an inverse relationship with silicic fraction except in Holocene. Below 15 ka BP, CaCO₃ dominates in sediments comprising more than 65%, such an increase is also seen in the coarse fraction. Coarse fraction from these sections contains abundant nodular type aggregates encrusting small forams. This period is marked by a high sedimentation rate comparable to Holocene. These parameters suggest that the productivity and precipitation have increased in the Holocene due to the intensification of the southwest monsoon. During the last glacial maximum and early deglacial period the high sedimentation rate indicates redeposition of the carbonates from the existing carbonate lithofacies situated between Quilon and Cape Comorin probably due to the slope instability.  相似文献